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Abstract 

In the world with axion the precise measurement of gravitational force and check of the 
equivalence principle at small distances, A ~ 1cm and less, could provide an additional test of 
CP symmetry. Using the chiral approach, we connect the strength of CP-odd axion-nucleon 
vertex with CP-odd pion-nucleon coupling for color electric dipole moments of quarks as a low 
energy source of CP violation. Strong limits on CP-odd coupling g n NN coming from atomic 
electric dipole moments show that gravitational experiments are six orders of magnitude 
away from the sensitivity level needed to detect axion-mediated forces. 
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1 Introduction 



The strong CP problem is one of the most intriguing issues of modern particle physics. The 
additional term in the QCD Lagrangian 

g 2 

violates P and CP symmetries ffl. In the electroweak theory, a diagonalization of quark 
mass matrices M u and Md involves chiral rotations and brings the additional contribution 
to theta term: 

6 = 6 + arg(detM u M d ) (2) 



Current experimental limits on the electric dipole moment (EDM) of the neutron put 
a severe constraint on 6 parameter. The chiral algebra calculation of the neutron EDM 
induced by the theta term || gives the following prediction: 

d n ~ 3.6 xl0- 16 6e- cm. (3) 

Together with the current neutron EDM constraints it implies the limit 6 < 3 • 10~ 10 . 

The puzzling smallness of 6 in comparison with a naive expectation 6 ~ 1 is usually 
referred to as strong CP problem. It could be solved theoretically in different manners. One 
can postulate that the starting point for ^-parameter is zero and mass matrices are hermitean 
at some high-energy scale due to some symmetry reasons [|J. The evaluation down to the 
electroweak scale must keep radiative corrections to 6 under the level of 10~ 9 . This scenario 
can hardly accommodate any sizable CP-violating phase beyond CP violation coming from 
Yukawa matrices (except, maybe, lepton sector). In all known models of CP violation 
which predict large EDMs and color EDMs (CEDMs) of quarks the radiative corrections 
to 6 give the EDM of the neutron already larger than current experimental bound. In 
the Minimal Supersymmetric Standard Model (MSSM), for example, this implies that a 
hidden dynamics responsible for supersymmetry breaking must respect CP to produce real 
soft-breaking parameters ||. 

In contrast, the existence of axion mechanism allows the dynamical relaxation of 6 || [7] 
and leaves a room for large CP violating effects induced by the effective operators dim> 5. 

The combination of CP-violation and PQ symmetry leads to a very interesting phe- 
nomenon, the long range interaction mediated by axion |§. It is shown schematically at Fig. 
1. The axion mass of 10 -5 eV corresponds to the characteristic length for this interaction 
of about 2 cm. It mimics gravity at these distances and therefore could be traceable in the 
measurement of gravitational forces at small distances ||. Moreover, this new interaction 
needs not respect equivalence principle and can show up in the precise Eotvos-Braginsky 
type of experiment in sub-cm region \W\. The latter is very promising because accuracy of 



this check at large distances can be as good as 10 



-12 
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In Standard Model, though, these effects are known to be negligibly small ||, |TTJ. In 
different modifications of supersymmetric models with new CP-violating sources besides the 
conventional Kobayashi-Maskawa phase there is a certain hope that the new long-range 
interaction can be detected by new precision measurements of gravitational forces at the 
small distances |12],[13|. Using the Naive Dimensional Analysis (NDA), the authors of work 



T3fl connected the strength of axion mediated force with maximal CP violation allowed by 
the limit on the electric dipole moment of neutron. Their result suffers from a rather large 
uncertainty (four orders of magnitude) and do not include the analysis of possible equivalence 
principle violation. In this letter we calculate the strength of the axion-matter coupling and 
dynamically induced theta term using chiral perturbation theory and some elements of the 
QCD sum rules. Our aim is to improve the accuracy of the predictions based on NDA and 
connect this vertex with CP-odd pion-nucleon coupling constant, strongly limited from the 
experiments searching for EDMs of diamagnetic atoms ]14| (For the detailed analysis of this 
coupling, see Refs. [IIR [UJ.) 



2 CP-odd interaction of the axion with matter 

Motivated by different supersymmetric models, we choose CEDMs as the most important 
effective operators and neglect possible effects from Weinberg three gluon operator and four- 
fermion type of interaction. Then the relevant part of the effective Lagrangian at the scale 
of 1 GeV is 

£e f/ = +~G%cT^ + i £ I^G^^ft + ... , (4) 

J a i=u,d,s 

were a is the axion field and f a is the axion decay constant. Dots stand here for the model 
dependent gradient interaction, d^a ■ j^, which can generate only contact interaction and 
therefore is irrelevant for our problem. 

In the absence of CP violation, nonremovable by PQ transformation, PQ symmetry sets 
theta parameter to zero @. The situation is different in the presence of extra CP- violating 
sources, communicated by the operators dim> 5. Thus, CEDMs of quarks drive theta 
parameter from zero to a value given by a ratio of two correlators: 

9 eff = where K = i \ f dxe ikx (0\T(^GG(x), ^GG(Q))\Q)\ (5) 

\K\ U 87T 57T ) k=0 

K 1 = i Udxe^(0\T(^GG(x), i Y,^9(Ga) l5 q l (0)\0) 



k=0 



and G a G a = GG, (Ga) = t a G a ^ u a^ u . The nonzero value of O e ff, in its turn, induces 
CP- violating axion- nucleon coupling ||: 

9aNN = °-f-{N\ ™* m ^ {Uu + dd + ss)\N). (6) 

f a m s m d + m s m u + m d m u 
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Even though we are already able to set a constraint on this vertex using (§), it is still 
instructive to calculate K\ and O e ff- 

The calculation of K is based on the use of the anomaly equation and the saturation of 
subsequent quark pseudoscalar correlator by light hadronic states JT7J. It yields a formula 
for the axion mass: 

m 2 = -— = m "f* mul7ld (7) 
a P a P a (m u + m d y { } 

Same technique can be exploited in the case of K\ (for the case of m u = m d the explicit 
derivation of K% can be found in Ref. Jl8[). Using the anomaly equation in the form proposed 
in 
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m d m s uj^ 5 u + m u m s d-f^ 5 d + m u m d s~f^ 5 s _ 

M m s m d + m s m u + m d m u 

2m u m d m s - _ a s ~ 

-{wy 5 u + d^d + S75SJ + — GG, 



m s m d + m s m u + m d m u " An 

we apply standard technique of current algebra. The correlator of interest, K±, can be 
rewritten in the form of the equal-time commutator, which we calculate easily, plus a term 
containing the singlet combination of pseudoscalars built from quark fields: 

K x = -i(0| (^ug{Ga)u+^dg{Ga)d + ^sg{Go)s\ |0> + (9) 

2 m s m d + m s m u + m d m u \m u m d m s I 



[ d 4 x(0\T{ im u m d m s ^ + ^ + s ^ x ^ ^ ^-qig(Gcr)n/5qi(0)}\0} 

J m s m d + m s m u + m d m u , 2 



The second term here is suppressed by an extra power of light quark masses in the numerator. 
It would bring a comparable contribution, though, if there were an intermediate hadronic 
state with a mass, vanishing in the chiral limit m, — > 0. At the same time the flavor structure 
of this term shows that the lightest intermediate state here is rf which is believed to remain 
heavy even if quark masses vanish. Thus the contribution from the second term is negligible 
in the limit m n m v r. 

Finally we get the dynamically induced theta term in the following compact form: 

e ^fi L+ d £+ iy (10) 

2 \m u m d m s J 

here is the ratio of the quark-gluon condensate to the quark condensate. It is known to 
sufficient accuracy from the baryon sum rules |2T 



2 (0\9q(Gcr)q\0) 2 
m n = — — — ~ U.8GeV . (11) 

Now we turn to the calculation of the axion-nucleon CP-violating coupling not related 
to 9. This coupling is given by the matrix element over the nucleon state of the structure 
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similar to K\\ 



ijdx{N\L ik *T(^GG(x),i £ f^ a (GV) 7 5*(0)l \N) (12) 

I 71 i=u,d,s J fc=Q 

Using the equation (^), we go along the same line as in the calculation of e // and connect 
Eq. (p!2|) with the matrix elements of the dipole-type operators over nucleon. The new thing 
here is the existence of additional contributions related with the diagrams of Fig. 2. Taking 
into account all contributions, after some algebra, we arrive at the following formula: 



m d m s m u 



m s m d + m s m u + m d m u 



_l(jV| [ ^Lu(Ga)u + —d(Ga)d + ^s(Ga)s) \N) + 
2 \m u m d m s I 



m 2 (N\huu - dd) du dd + (uu + dd - 2ss) du + dd 2 f s \N) 
UN '2 V J m u + m d K ' m u + m d + 4m s 1 ' 



(13) 



It should be noted here that the CP-odd nature of the interaction aNN was crucial for the 
transformation of Eq. (|12"D to the series of matrix elements from local operators (0). It 
allows to drop all "nonlocal" terms proportional to fc Ai (A^|T(g7 At 7 5 g(x), iqg{Ga)^^q{Q)\N) as 
leading to k 2 in the numerator. It would be not the case, for example, if we were calculating 
corrections to CP-even interaction iaN , y 5 N caused by qg(Ga)q where similar matrix elements 
are not negligible. 



Further progress is related with the application of low energy theorem [21 in + channel 
to the matrix element of the quark-gluon operator 

(N\qg(Ga)q\N) ~ j^(jV|gg|iV>. (14) 



For (N\qq\N) we have the following numbers for certain flavors |2J 

(p\uu\p) ~ 4.8; (p\dd\p) ~ 4.1; (p\ss\p) ~ 2.8 (15) 

with the obvious generalization on the neutron. 

From these relations we are able to learn one unfortunate thing. The dynamics of strong 
interaction makes CP-violating coupling of axion with neutron and proton be nearly equal 
even if CP violation at small distances is isospin asymmetric (exept, maybe, some very 
specific cases), 

9app 9ann ^ 9app + 9 arm = ^9 aNN (16) 



The strength of isospin-symmetric coupling g a NN follows from Eq. ([13|) after the substitution 
of©-(0), 

1 m u m d 2 / d u d d n -ds\ n -x 

g a NN ^ - f ; m 3.8 + — - 0.7 — 17 

f a m u + m d \ m u m d m s J 

At this point we should plug in the limits on di coming from the EDM-type of experiments. 
As it was shown in [|l|, [16[ , the data on atomic EDMs set better limits on d u and d d than 
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neutron EDM does. For us, however, it is more important that g a NN can be connected with 
pion-nucleon CP-odd vertex with little degree of uncertainty. This vertex, g n NN, leads to the 
T-odd nucleon-nucleon interaction NNN'i^N' , to Schiff moment of nuclei and ultimately to 



an atomic EDM. This calculation was done in Refs. [[L5|, [16[ by means of the same technique 
that we exploit here and we simply quote this result: 

! — ~ — 1 ) ~ — 1111 ~\~ cfn 

QttNN = —(N\d u ug(Ga)u - d d dg(Ga)d\N) ~ --rn^K - d d )(N\^—\ N). (18) 
2/tt U 6 2 

Again, (N\uu — dd\N) was neglected in comparison with (N\uu + dd\N). Formulae (|17|)- 
(|18|) solve, in principle, the problem of correspondence between g^NN and g a NN- However, 
instead of writing a general formula connecting these two quantities, it is convenient to go to 
some specific and more instructive physical cases making some assumptions about underlying 
CP-violating physics. 

Case 1. Proportionality 
By "proportionality" we understand a certain situation when d u /m u = d d /m d = d s /m s . 
It happens, for example, in MSSM when the relative CP-violating phase between trilinear 
soft-breaking parameter and gaugino mass is mediated by gluino exchange diagram. 

In this case the relation between two couplings is given by the following formula: 

f n 2m u m d 
fa mj - mi 

As to the potential violation of the equivalence principle by this interaction, it is given by 
the following relation: 

gann ~ fi'app _ (p\uU - dd\p) 



gann + gapp (p\uU + dd\p) 



0.08. (20) 



Case 2. One-flavor dominance 
This is an opposite to the previous case situation when di/mi is enhanced for certain flavor. 
This may happen within different Grand Unification scenarios 



For u(d)-dominance the connection between g a NN and g n NN is unambiguous; it follows 
from general formulae flTB| ) and ( JTB| ) before taking matrix elements, so that major QCD- 
related uncertainties are gone. The resulting relation takes the simplest possible form: 

U m d(u) / 01 x 

g a NN — g-KNN-^ ; , l^-LJ 

f a m u + m d 

where the subscripts d(u) correspond to the case of up(down)-flavor dominance. 

If d u —d d is close to zero (although, we do not know explicit examples of such a model), the 
coupling g-xNN becomes small. It does not mean, of course, that g a NN can grow uncontrollably 
on account of the limit from neutron EDM which is sensitive to d u + d d combination ||16||. 
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3 Comparison with experimental data 



The experimental bound on EDM of 199 Hg (|TJJ translates into the following stringent limit 
on g-xNN 

Qirpp ^ 9nNN < 2 ■ KT 11 (22) 

In the most interesting case of proportionality we get the following bound on the ratio 
of axion-mediated long-range interaction to gravitational force for the range of new force 
A = m~ l : 

^™ = = 5 " lO 11 ^ (^P) 2 < 2 • 10- 10 (^) 2 . (23) 

Fgravity 4:71^171% \ X J \ A / 

This ratio is tremendously small as compared to the current experimental sensitivity at the 
level of lO" 4 . 

More promising in terms of experimental accuracy is a high-precision check of the equiv- 
alence principle at small distances. The existing experiment sets the limit on the relative 
accelerations of Cu and Pb toward a uranium attractor as precise as 10~ 8 at the distances 
A > 10 cm |fTof| . What can we expect from the interaction mediated by axion in this case? 
For the case of maximal isospin-asymmetric CP-violation (<i u -dominance) the estimate of 
the equivalence principle violation by axion mediated forces takes the following form: 



acu - apb ^ (N - Z) Pb (p\uu - dd\p) F axion ^ o in _ 12 (lcm\ 

dgr-avity 2A Pb (p\lLU + dd\p) F gravity 

The results (EBP- (El) are summarized in Fig 3. 
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< 2 • KT 12 f — J . (24) 



4 Conclusions 



In the case of CP violation communicated by color EDMs, the axion-nucleon CP-violating 
vertex is calculated by means of the current algebra. We have shown that the sensitivity 
of current gravitational experiments at small, ~ 1cm, distances needs to be improved by 
six orders of magnitude to probe the CP-violating interaction of axion with matter. Our 
conclusion is very firm and do not contain a large degree of uncertainty reported in the 
previous analysis [[L3|]. It is based on the relations between g a NN and g^NN, Eqs. (|19"D- 



(pl|), which are held to better accuracy than the calculation of m a , g a NN and g n NN itself. 
Although, color EDMs are, perhaps, the most interesting CP-violating operators, one can 
implement similar procedure for all variety of CP-odd operators dim=6 built from the quark 
fields. For Weinberg operator as a source, though, we are not able to give a better estimate 
than that based on NDA. 

Our result for F axion / F grav i ty , Eq. fl23|) , is significantly lower than the central value of the 
prediction based on NDA and neutron EDM data ||13||. Once again it shows usefulness of 



the limits on CP- violation extracted from EDM of i99 Hg [14 
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Figure captions. 

Fig. 1 The symbolic picture of axion exchange between two nucleons with CP violation 
in both vertices. 

Fig. 2 Additional contribution to g a NN, not related to 6 or (N\qg(Ga)q\N) . 

Fig. 3a The comparison of maximal strength for F axion / F gravity (straight line) with current 
experimental sensitivity (curve in the upper corner) ||]. 

Fig. 3b The comparison of predicted maximal strength for Aa/a with the limits from 
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Figure 2. 



Log(F/F) 
a g 



Log(A a/a gr ) 





Figure 3 a. 



Figure 3b. 
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